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We consider an electrostatic qubit, interacting with fluctuating charge of a single electron tran- 
sistor (SET) in the framework of an exactly solvable model. The SET plays role of an environment 
affecting the qubits' parameters in a controllable way. We derive the rate equations describing the 
dynamics of the entire system for an arbitrary qubit-SET coupling. Solving these equations we 
obtain decoherence and relaxation rates of the qubit, as well as the spectral density of qubit param- 
eters' fluctuations. We found that in a weak coupling regime decoherence and relaxation rates are 
directly related to the spectral density taken at either zero or Rabi frequency, depending on which 
qubit parameter is fluctuating. In the latter case our result coincides with that of the spin-boson 
model in the weak coupling limit, despite different origin of the fluctuations. We show that this 
relation holds also in the presence of weak back-action of the qubit on the environment. In case of 
strong back-action such a simple relationship no longer holds, even if qubit-SET coupling is small. 
It does not hold also in the strong coupling regime, even in the absence of the back-action. In 
addition, we found that our model predicts localization of the qubit in the strong-coupling regime, 
resembling that in the spin-boson model. 
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I. INTRODUCTION 

The influence of environment on a single quantum sys- 
tem is the issue of crucial importance in quantum in- 
formation science. It is mainly associated with decoher- 
ence, or dephasing, which transforms any pure state of 
a quantum system into a statistical mixture. Despite a 
large body of theoretical work devoted to decoherence, its 
mechanism has not been clarified enough. For instance, 
how decoherence is related to environmental noise, in 
particular in the presence of back-action of the system 
on the environment (quantum measurements). More- 
over, decoherence is often intermixed with relaxation. Al- 
though each of them represents an irreversible process, 
decoherence and relaxation affect quantum systems in 
quite different ways. 

In order to establish a relation between the fluctua- 
tion spectrum and decoherence and relaxation rates one 
needs a model that describes the effects of decoherence 
and relaxation in a consistent quantum mechanical way. 
An obvious candidate is the spin-boson model^'^ which 
represents the environment as a bath of harmonic oscil- 
lators at equilibrium, where the fluctuations obey Gaus- 
sian statistics^. Despite its apparent simplicity, the spin- 
boson model cannot be solved exactly^. Also, it is hard 
to manipulate the fluctuation spectrum in the framework 
of this model. In addition, mesoscopic structures may 
couple only to a few isolated fluctuators, like spins, lo- 
cal currents, background charge fluctuations, etc. This 
would require models of the environment, different from 
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FIG. 1: Electrostatic qubit, realized by an electron trapped 
in a coupled-dot system (a) , and its schematic representation 
by a double- well (b). ilo denotes the coupling between the 
two dots. 

the spin-boson model (see for instace^'^i^ii^'^'iSiiii'i^) . In 

general, the environment can be out of equilibrium, like 
a steady-state fluctuating current, interacting with the 
qubilji^iii^'i^ii^. This for instance, takes place in the con- 
tinuous measurement (monitoring) of quantum systems^i 
and in the "control dephasing" experiment o^^i^^'^" . All 
these types on non-Gaussian and non-equilibrium envi- 
ronments attracted recently a great deal of attention^^. 

In this paper we consider an electrostatic qubit, which 
can be viewed as a generic example of two-state systems. 
It is realized by an electron trapped in coupled quantum 
dots^^i^i^. Fig. 1. Here Ei and E2 denote energies of 
the electron states in each of the dots and Hq is a cou- 
pling between these states. It is reasonable to assume 
that the decoherence of a qubit is associated with fluctu- 
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FIG. 2: Qubit near Single Electron Transistor. Here i5;,r and 
Eo denote the energy levels in the left (right) reservoirs and in 
the quantum dot, respectively, and fiL,R are the correspond- 
ing chemical potentials. The electric current I generates fluc- 
tuations of the electrostatic opening between two dots (a), or 
it fluctuates the energy level of the nearest dot (b). 



ations of the qubit parameters, Ei^2 and flo, generated 
by the environment. Indeed, a stochastic averaging of 
the Schrodinger equation over these fluctuations param- 
eters results in the qubit's decoherence, which transfers 
any qubit state into a statistical mixtur o^^i^^ . In general, 
one can expect that the fluctuating environment should 
result in the qubit's relaxation, as well, as for instance in 
the phenomenological Redfield's description of relaxation 
in the magnetic resonance^^. 

As a quantum mechanical model of the environment 
we consider a Single Electron Transistor (SET) capac- 
itively coupled to the qubit, e.g., Fig. 2. Such setup 
has been contemplated in numerous solid state quantum 
computing architectures where SET plays role of a read- 
out devicei^iiiiSS,!^ and contains most of the generic fea- 
tures of a fluctuating non-equilibrium environment. The 
discreteness of the electron charge creates fluctuations 
in the electrostatic field near the SET. If the electro- 
static qubit is placed near the SET, this fluctuating field 
should affect the qubit behavior as shown in Fig. 2. It 
can produce fluctuations of the tunneling coupling be- 
tween the dots (off-diagonal coupling) by narrowing the 
electrostatic opening connecting these dots, as in Fig. 2a, 
or make the energy levels of the dots fluctuate, as shown 
schematically in Fig. 2b. Note that while in some regimes 
the SET operates as a measuring devicei^iii, in other 
regimes it corresponds purely to a source of noise. In- 
deed, if the energy level Eq, Fig. 2, is deeply inside 
the voltage bias - the case we consider in the begin- 
ning, the SET current is not modulated by the qubit 
electron. In this case the SET represents only the fluc- 
tuating environment affecting the qubit behavior ("pure 
environment"—) . 

A similar model of the fluctuating environment has 
been studied mostly for small bias (linear response) or 



for the environment in an equilibrium. Here, however, 
we consider strongly non-equilibrium case where the bias 
voltage applied on the SET (y = /i^ — /i/;) is much 
larger than the levels widths and the coupling between 
the SET and the qubit. In this limit our model can be 
solved exactly for both weak and strong coupling (but 
is still smaller than the bias voltage). This constitutes 
an essential advantage with regard to perturbative treat- 
ments of similar models. For instance, the results of our 
model can be compared in different regimes with phe- 
nomenological descriptions used in the literature. Such 
a comparison would allow us to determine the regions 
where these phenomenological models are valid. 

Since our model is very simple in treatment, the deco- 
herence and relaxation rates can be extracted from the 
exact solution analytically, as well as the time-correlator 
of the electric charge inside the SET. This would make 
it possible to establish a relation between the frequency- 
dependent fluctuation spectrum of the environment and 
the decoherence and the relaxation rates of the qubit, 
and to determine how far this relation can be extended. 
We expect that such a relation should not depend on 
a source of fluctuations. This point can be verified by 
a comparison with a similar results obtained for equi- 
librium environment in the framework of the spin-boson 
modelii^. 

It is also important to understand how the decoher- 
ence and relaxation rates depend on the frequency of 
the environmental fluctuations. This problem has been 
investigated in many phenomenological approaches for 
"classical" environments at equilibrium. Yet, there still 
exists an ambiguity in the literature related to this point 
for non-equilibrium environment. For instance, it was 
found by Levinson that the decoherence rate, generated 
by fluctuations of the energy level in a single quantum dot 
is proportional to the spectral density of fluctuations at 
zero frequencj«21. The same result, but for a double-dot 
system has been obtained by Rabenstein et al?^ . On 
the other hand, it follows from the Redfied's approach 
that the corresponding decoherence rate is proportional 
to the spectral density at the frequency of the qubit's 
oscillations (the Rabi frequency)^. Since our model is 
the exactly solvable one, we can resolve this ambiguity 
and establish the appropriate physical conditions that 
can result in different relations of decoherence rate to 
the environmental fluctuations. 

The most important results of our study are related to 
the situation when back-action of the qubit on the envi- 
ronment takes place. This problem did not receive such 
a considerable amount of attention in the literature as, 
for example, the case of "inert" environment. This is in 
spite of a fact that the back-action always takes place in 
the presence of measurement. There are many questions 
related to the effects of a back-action. For instance, what 
would be a relation between decoherence (relaxation) of 
the qubit and the noise spectrum of the environment? 
Or, how decoherence is affected by a strong response of 
environment? We believe that our model appears to be 
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more suitable for studying these and other problems re- 
lated to the back-action than most of the other existing 
approaches. 

The plan of this paper is as follows: Sect. II presents a 
phenomenological description of decoherence and relax- 
ation in the framework of Bloch equations, applied to the 
electrostatic qubit. Sect. Ill contains description of the 
model and the quantum rate-equation formalism, used 
for its solution. Detailed quantum-mechanical derivation 
of these equations for a specific example is presented in 
Appendix A. Sect. IV deals with a configuration where 
the SET can generate only decoherence of the qubit. We 
consider separately the situations when SET produces 
fluctuations of the tunneling coupling (Rabi frequency) 
or of the energy levels. The results are compared with 
the SET fluctuation spectrum, evaluated in Appendix B. 
Sect. V deals with a conflguration where the SET gener- 
ates both decoherence and relaxation of the qubit. Sect. 
VI is summary. 



II. DECOHERENCE AND RELAXATION OF A 
QUBIT 

In this section we describe in a general phenomeno- 
logical framework the effect of decoherence and relax- 
ation on the qubit behavior. Although the results are 
known, there still exists some confusion in the literature 
in this issue. We therefore need to deflne precisely these 
quantities and demonstrate how the corresponding deco- 
herenece and relaxation rates can be extracted from the 
qubit density matrix. 

Let us consider an electrostatic qubit, realized by an 
electron trapped in coupled quantum dots. Fig. 1. This 
system is described by the following tunneling Hamilto- 
nian 

Hqh = Eia\ai + E2a\a2 ~ VLQ{a\ai + a\a2) (1) 

where a\ 2,01.2 are the creation and annihilation opera- 
tors of the electron in the flrst or in the second dot. For 
simplicity we consider electrons as spinless fermions. In 
addition, we assume that a\ai -f a\a2 = 1, so that only 
one electron is present in the double-dot. The electron 
wave function can be written as 



!*(<))= [6«(i)al+ 6(2) (i)4 
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(2) 



where ^(^'^^(i) are the probability amplitudes for find- 
ing the electron in the first or second well, obtained 
from the Schrodinger equation i9t| 5* (f)) — Hqh\^{t)) (we 
adopt the units where h = \ and the electron charge 
e — 1). The corresponding density matrix, (Jjji{t) = 
b'^i){t)b^'y)*{t), with = {1,2}, is obtained from the 
equation idt a = [i/ , cr] . This can be written explicitly as 

0-11 = irJo(a'2i - a'12) (3a) 
(j\2 = -iecri2 + irio(l - 2crii) , (3b) 



where 0-22 (t) = 1 - crii(i), (J2\{t) = cr j['2(0 and e = 
E\ — £2- Solving these equations one easily finds that 
the electron oscillates between the two dots (Rabi oscil- 
lations) with frequency lub. = ■\/4ilQ + e^. For instance, 
for the initial conditions crii(O) = 1 and (Ti2(0) — 1, the 
probability of finding the electron in the second dot is 
C22(i) = 2(ilo/^7?)'^(l ~ cosujjit). This result shows that 
for e ^ rig the amplitude of the Rabi oscillations is small, 
so the electron remains localized in its initial state. 

The situation is different when the qubit interacts with 
the environment. In this case the (reduced) density ma- 
trix of the qubit a{t) is obtained by tracing out the en- 
vironment variables from the total density matrix. The 
question is how to modify Eqs. ([3]), written for an isolated 
qubit, in order to obtain the reduced density matrix of 
the qubit, cr(t). In general one expects that the environ- 
ment could affect the qubit in two different ways. First, it 
can destroy the off-diagonal elements of the qubit density 
matrix. This process is usually referred to as decoherence 
(or dephasing). It can be accounted for phcnomenolog- 
ically by introducing an additional (damping) term in 
Eq. pb)l . 



-0'12 



(4) 



where is the decoherence rate. As a result the qubit 
density- matrix ait) becomes a statistical mixture in the 
stationary limit. 



1/2 
1/2 



(5) 



This happens for any initial conditions and even for large 
level displacement, e ^ iloiTd (provided that JIq 7^ 0). 
Note that the statistical mixture ([5]) is proportional to 
the unity matrix and therefore it remains the same in 
any basis. 

Secondly, the environment can put the qubit in its 
ground state, for instance via photon or phonon emis- 
sion. This process is usually referred to as relaxation. 
For a symmetric qubit we would have 



ait) 



1/2 1/2 
1/2 1/2 



(6) 



In contrast with decoherence, Eq. ([5]), the relaxation pro- 
cess puts the qubit into a pure state. That implies that 
the corresponding density matrix can be always written 
as SiiSij in a certain basis (the basis of the qubit eigen- 
states). This is in fact the essential difference between 
decoherence and relaxation. With respect to elimination 
of the off-diagonal density matrix elements, note that re- 
laxation would eliminate these terms only in the qubit 's 
eigenstates basis. In contrast, decoherence eliminates the 
off-diagonal density matrix element in any basis (Eq. 

In fact, if the environment has some energy, it can put 
the qubit into an exited state. However, if the qubit is 
finally in a pure state, such excitation process generated 
by the environment affects the qubit in the same way as 
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relaxation: it eliminates the ofF-diagonal density matrix 
elements only in a certain qubit's basis. Therefore exci- 
tation of the qubit can be described phenomenologically 
on the same footing as relaxation. 

It is often claimed that decoherence is associated with 
an absence of energy transfer between the system and 
the environment, in contrast with relaxation (excitation). 
This distinction is not generally valid. For instance, if the 
initial qubit state corresponds to the electron in the state 
\E2), Fig. 1, the final state after decoherence corresponds 
to an equal distribution between the two dots, (E) — 
{El + i?2)/2. In the case of Ei » this process would 
require a large energy transfer between the qubit and the 
environment. Therefore decoherence can be consistently 
defines as a process leading to a statistical mixture, where 
all states of the system have equal probabilities (as in 
Eq. ©). 

The relaxation (excitation) process can be described 
most simply by diagonalizing the qubit Hamiltonian, 
Eqs. HI), to obtain i/qb = £'+0^0+ + £'_a^a_, where 



the operators a± are obtained by the corresponding ro- 
tation of the operators ai ^SS,. Here E+ and E_ are the 
ground (symmetric) and excited (antisymmetric) state 
energies. Then the relaxation process can be described 
phenomenologically in the new qubit basis |±) = a]j_|0) 
as 

&—{t) = -Tr(7—{t) (7a) 
&+-it) = i{E- - E+)a+-{t) - Y^+-(t) , (7b) 

where cr_|_+(t) = 1 — (T__(t), a ^^{t) — (T^„(t) and F^ is 

the relaxation rate. 

In order to add decoherence, we return to the orig- 
inal qubit basis |1,2) = a|2|0) and add the damping 
term to the equation for the off-diagonal matrix elements, 
Eq. (HI) . We arrive at the quantum rate equation describ- 
ing the qubit's behavior in the presence of both decoher- 
ence and relaxation^SiS 



&11 — ifio(o'2i — (T12) (o'i2 -I- (J21) ^ 



(2aii-l) + r,- 



CT12 — ~ifJi2 + 



■0 , T- '^^ 



(l-2aii) + F, 



K - 2 '^12 - «^(a-i2 + 0-21) 



y^i2 



(8a) 
(8b) 



where e = (e^ + 4rig)^/^ and k ~ Vlo/l. In fact, these 
equations can be derived in the framework of a particu- 
lar model, representing an electrostatic qubit interacting 
with the point-contact detector and the environment, de- 
scribed by the Lee model HamiltonianS^. 

Equations (HI) can be rewritten in a simpler form by 
mapping the qubit density matrix a = {(J 11, cr 12, 21} to 
a "polarization" vector S{t) via a{t) = [1 + t ■ S{t)]/2, 
where r^^y.z are the Pauli matrices. For instance, one 
obtains for the symmetric case, e = 0, 



Sz — — 2'^^ ^ '^^^ 
Sy — 2^1q Sz Srii 



■'y — 2 
Fd + 2F 



St. — 



{,Sx Sx) 



(9a) 
(9b) 
(9c) 



where Sx — Sx{t ^ 00) = 2Vr/ (F<i-l-2Fr). One finds that 
Eqs. ^ have a form of the Bloch equations for spin- 
precession in the magnetic fieldS^, where the effect of 
environment is accounted for by two relaxation times for 
the different spin components: the longitudinal Ti and 
the transverse T2, related to F^; and 2F,. as 



and T2^ = 



Td+Tr 



(10) 



The corresponding damping rates, the so-called "depolar- 
ization" (Fi = 1/Ti) and the "dephasing" (F2 = I/T2) 
are used for phenomcnological description of two-level 
systems^. However, neither Fi nor F2 taken alone would 
drive the qubit density matrix into a statistical mixture 
Eq. (O or into a pure state Eq. 

In contrast, our definition of decoherence and relax- 
ation (excitation) is associated with two opposite effects 
of the environment on the qubit: the first drives it into 
a statistical mixture, whereas the second drives it into 
a pure state. We expect therefore that such a natural 
distinction between decoherence and relaxation would be 
more useful for finding a relation between these quantities 
and the environmental behavior than other alternative 
definitions of these quantities existing in the literature. 

In general, the two rates, F^ ,., introduced in phe- 
nomcnological equations (I5|) , (O , are consistent with our 
definitions of decoherence and relaxation. The only ex- 
ception is the case of F^ = and = 0, where are no 
transitions between the qubit's states even in the pres- 
ence of the environment ("static" qubit). One easily finds 
from Eqs. (|3a|) . ([4|) that ai2{t) ^ for i ^ 00, whereas 
the diagonal density- matrix elements of the qubit remain 
unchanged (so-called "pure dephasing"i^i^) : 



a(t) 



f7ll(0) 
CT22(0) 



(11) 
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Thus, if the initial probabilities of finding the qubit in 
each of its states are not equal, crii(O) ^ 0-22(0), then the 
final qubit state is neither a mixture nor a pure state, but 
a combination of the both. It implies that Yd in Eqs. ([5]) 
would also generate relaxation (excitation) of the qubit. 
Note that in this case the off-diagonal density-matrix el- 
ements, absent in Eq. ljlip . would reappear in a different 
basis. This implies that the "pure dephasing" ^'^"^ occurs 
only in a particular basis. 

Let us evaluate the probability of finding the electron 
in the first dot, (J\\{t). Solving Eqs. ^ for the initial 
conditions 0-11(0) = 1, 012(0) = 0, we find^: 

^ii(t) - ^ + '—^ (C^ie-^-* + ^26-^+*) (12) 

where e± = i(rd ± f^), = - 64rj^ and Ci,2 = 
l±(rd/f2). Solving the same equations in the limit of i — > 
00, we find that the steady-state qubit density matrix is 

nif\'=:^i rj(r, + 2r,)\ 

''^^^ \ r,./(r<i + 2r,) 1/2 ) ■ ^^"^^ 

Thus the off-diagonal elements of the density matrix 
can provide us with a ratio of relaxation to decoherence 
rates^. 



III. DESCRIPTION OF THE MODEL 

Consider the setup shown in Fig. 2. The entire system 
can be described by the following tunneling Hamiltonian, 
represented by a sum of the qubit and SET Hamiltonians 
and the interaction term, H = Hqh + -ffsET + ^int- Here 
ffqb is given by Eq. ^ and describes the qubit. The sec- 
ond term, Hset, describes the single-electron transistor. 
It can be written as 

Hset = ^ Eicjci + ^ Erclcr + EqcIcq 

I r 

+ ^{nic\co + nrclco + H.C.) , (14) 



I 



where h^^\t), ba\t) are the probability amplitudes to 
find the entire system in the state described by the cor- 



where cj ^ and Ci^r are the creation and annihilation elec- 
tron operators in the state Ei^r of the right or left reser- 
voir; Cq and Co are those for the level Eg inside the quan- 
tum dot; and Qi r are the couplings between the level Eq 
and the level Ei^r in the left (right) reservoir. In order to 
avoid too lengthy formulaes, our summation indices I, r 
indicate simultaneously the left and the right leads of the 
SET, where the corresponding summation is carried out. 
As follows from the Hamiltonian (fT4|) . the quantum dot 
of the SET contains only one level (Eq). This assumption 
has been implied only for the sake of simplicity for our 
presentation, although our approach is well suited for a 
case of n levels inside the SET, SqcJco J2n -^'"■'^I'^ri, 
and even when the interaction between these levels is in- 
cluded (providing that the latter is much less or much 
larger than the bias V)^^^. We also assumed a weak 
energy dependence of the couplings fli^r — ^l,r- 

The interaction between the qubit and the SET, Hint, 
depends on a position of the SET with respect to the 
qubit. If the SET is placed near the middle of the qubit. 
Fig. 2a, then the tunneling coupling between two dots of 
the qubit in Eq. ([1]) decreases, flo — SQq, whenever 

the quantum dot of the SET is occupied by an electron. 
This is due to the electron's repulsive field. In this case 
the interaction term can be written as 

Hint = Snclco{a\a2 + alai) . (15) 

On the other hand, in the configuration shown in Fig. 2b 
where the SET is placed near one of the dots of the qubit, 
the electron repulsive field displaces the qubit energy lev- 
els by AE — U. The interaction terms in this case can 
be written as 

Hint = U alaiclcQ . (16) 

Consider the initial state where all the levels in the left 
and the right reservoirs are filled with electrons up to the 
Fermi levels hl.r respectively. This state will be called 
the "vacuum" state |0). The wave function for the entire 
system can be written as 



(17) 

I 

responding creation and annihilation operators. These 
amplitudes are obtained from the Schrodinger equation 



6(1) it)a\ + ^ fc^}) {t)a\clci + ^ 6« it)a\ 



Lr 



clci+ b''oJn,{t)a\clclciCi> 

l<l',r 



-6(^)(i)4 + ^6(f(t)4cSQ-f ^6(f(i)4ctQ+ J2 ^owrW4c^^ 



Ci-ClCl' 



Lr 



l<l',r 
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i\'^{t)) = H\^{t)), supplemented with the initial condi- 
tion 6(i)(0) = pi, 6(2)(o) = and b^J\o) = 0, where 
Pi_2 are the amplitudes of the initial qubit state. 

Note that Eq. (fT7| implies a fixed electron number {N) 
in the reservoirs. At the first sight it would lead to deple- 
tion of the left reservoir of electrons over the time. Yet 
in the limit of A'^ ^ oo (infinite reservoirs) the dynamics 
of an entire system reaches its steady state before such a 
depletion takes place^i^^. 



0-aai.t) 
<Jdd{t) 



In was shown in^ii^ that the trace over the reservoir 
states in the system's density matrix can be performed 
in the large bias limit (strong non-equilibrium limit) 

T^ = A«L-Mi?»r,r!o,C/ (19) 

where the level (levels) of the SET carrying the current 
are far away from the chemical potentials, and F is the 
width of the level Eq- In this derivation we assumed 
only weak energy dependence of the transition ampli- 
tudes fli^r = ^L,Ft, and the density of the reservoir states, 
p(Ei^r) ~ pL.R- As a result we arrive at Bloch-type rate 
equations for the reduced density matrix without any 
additional assumptions. The general form of these equa- 
tions is^i2^ 

ajf = i{Ey - Ej)ajj' + «X1 {'^jk^k^r - f^j^fco-fej') 

k 

— '^2^V2T^ p{(T'jk^k~,k'^k'^]' + (Tk]'^k~,k'^k'^j) 
k,k' 

+ ^ 'P27rp {ilk^jftk'^]' + flk^j'^ik'^j)akk' (20) 

k,k' 

Here rifc_,/c' denotes the single-electron hopping ampli- 
tude that generates the k ^ k' transition. We distinguish 
between the amplitudes 57 describing single-electron hop- 
ping between isolated states and describing transitions 
between isolated and continuum states. The latter can 
generate transitions between the isolated states of the 
system, but only indirectly, via two consecutive jumps of 
an electron, into and out of the continuum reservoir states 
(with the density of states p) . These transitions are rep- 
resented by the third and the fourth terms of Eq. 
The third term describes the transitions {k ^ k' ^ j) 
or (k —f k' j'), which cannot change the number of 



The behavior of the qubit and the SET is given by 
the reduced density matrix, ass' (t) ■ It is obtained from 
the entire system's density matrix |\E'(t)) (\['(i)| by trac- 
ing out the (continuum) reservoir states. The space of 
such a reduced density matrix consists of four discrete 
states s, s' — a, 6, c, d, shown schematically in Fig. 3 for 
the setup of Fig. 2a. The corresponding density-matrix 
elements are directly related to the amplitudes b{t), for 
instance. 



(18a) 
(18b) 
(18c) 

I 

electrons in the collector. The fourth term describes the 
transitions (fc — s- j and k' — > j') or (fc — s- j' and k' — > j) 
which increase the number of electrons in the collector 
by one. These two terms of Eq. (|^D|) are analogues of the 
"loss" (negative) and the "gain" (positive) terms in the 
classical rate equations, respectively. The factor 7^2 — ±1 
in front of these terms is due anti-commutation of the 
fermions, so that 7^2 = —1 whenever the loss or the gain 
terms in Eq. (|20p proceed through a two-fermion state of 
the dot. Otherwise 7^2 = 1- 

Note that the reduction of the time-dependent 
Schrodinger equation, z|^'(t)) = H\'^{t)), to Eqs. ^ 
is performed in the limit of large bias without explicit 
use of any Markov-type or weak coupling approxima- 
tions. The accuracy of these equations is respectively 
max(r, Qq, U, T)/\pl,r ~ Ej\. A detailed example of this 
derivation is presented in Appendix A for the case of res- 
onant tunneling through a single level. The derivation 
there and in Refsi^i^^ were performed by assuming zero 
temperature in the leads, T = 0. Yet, this assumption 
is not important in the case of large bias, providing the 
levels carrying the current are far away from the Fermi 
energies, \pl,r - > T. 



IV. NO BACK-ACTION ON THE 
ENVIRONMENT 

A. Fluctuation of the tunneling coupling 

Now we apply Eqs. (^0]) to investigate the qubit's be- 
havior in the configurations shown in Fig. 2. First we 
consider the SET placed near the middle of the qubit. 
Figs. 2a, 3. In this case the electron current through the 
SET will influence the coupling between two dots of the 
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(a) 



LRi^ LTH^ LTHL 



(b) 



(c) 



(d) 



FIG. 3: The available discrete states of the entire system cor- 
responding to the setup of Fig. 2a. Vl,r denote the tunneling 
rates to the corresponding reservoirs and fig = f^o — 5^- 



qubit, making it fluctuate between the values f2o a-nd 
f2g = r^o — 5^1. The corresponding rate equations can be 
written straightforwardly from Eqs. ([20]) . One finds, 



CTaa = 










(21a) 


CTbf) = 


-rijCTfcb + 


^L(^aa - 


in'oiabd 


— CTdb), 


(21b) 


CTcc = 


-Tidcc + 


^R(7dd — 






(21c) 


O'dd = 




- TlCJcc — 




- 0-bd), 


(21d) 


CTac = 


— ieoo'ac — 




- <^cc) - 












H 




(21e) 


O'bd = 


— ^eo^bd — 


iQQ{abb - 














H 


- riCTac, 


(21f) 



where T^ n — 2Tr\nL j^\^pL fi are the tunneling rates from 
the reservoirs and eg = Ei — E2- 

These equations display explicitly the time evolution 
of the SET and the qubit. The evolution of the for- 
mer is driven by the first two terms in Eqs. (l2Ta|) - (pTdl) . 
They generate charge-fluctuations inside the quantum 
dot of the SET (the transitions a< — >b and c< — >d), 
described by the "classical" Boltzmann-type dynamics. 
The qubit 's evolution is described by the Bloch-type 
terms (c.f. Eqs. ([2])), generating the qubit transitions 
(a< — >c and b< — >d). Thus Eqs. ((2T|) are quite general, 
since they described fluctuations of the tunneling cou- 
pling driven by the Boltzmann-type dynamics. 

The resulting time evolution of the qubit is given by 
the qubit (reduced) density matrix: 



O'll(i) = (Taa{t) + O-bb(t) , 
O'l2(0 = <yac{t) + (Tbd{t) , 



(22a) 
(22b) 



and (T22(i) = 1 - o■ll(^)■ 
Similarly, the charge fluctuations of SET are deter- 
mined by the probability of finding the SET occupied. 



Piit) = cTbb{t) + add{t) . 
It is given by the equation 

A(i) = ri-rPi(t), 



(23) 



(24) 



obtained straightforwardly from Eqs. (pij) . Here F = 
F^ -I- Ffl is the total width. The same equation for Pi{t) 
can be obtained if the qubit is decoupled from the SET 



{50. = 0). Thus there is no back-action of the qubit on 
the charge fluctuations inside the SET in the limit of 
large bias voltage. 

Consider first the stationary limit, t oo, where 
Pi{t) and &{t) 0. It follows from Eq. ^ that the 
probability of finding the SET occupied in this limit is 
Pi = Fi/F. This implies that the fluctuations of the cou- 
pling r^Oi induced by the SET, would take place around 
the average value fl = Hq — Pi Sfl. 

With respect to the qubit in the stationary limit, one 
easily obtains from Eqs. pT|) that the qubit density ma- 
trix always becomes the statistical mixture ([5]), when 
t oo. This takes place for any initial conditions and 
any values of the qubit and the SET parameters. There- 
fore the effect of the fluctuating charge inside the SET 
does not lead to relaxation of the qubit, but rather to its 
decoherence. 

It is important to note, however, that for the aligned 
qubit, e = 0, the decoherence due to fluctuations of the 
tunneling coupling Oq is not complete. Indeed, it follows 
from Eqs. that d/dt[Re ai2{t)] = 0. The reason is 
that the corresponding operator, a\a2 + a\ai commutes 
with the total Hamiltonian H = Hq^ + Hset + ^^int, 
Eqs. HD, dHD and (15), for Ei = E2. As a result, 
Recri2(i) = Re (712(0). 

In order to determine the decoherence rate analytically, 
we perform a Laplace transform on the density matrix, 
a{E) = a{t) exp{-iEt)dE. Then solving Eq. ^ we 
can determine the decoherence rate from the locations 
of the poles of a{E) in the complex £'-plane. Consider 
for instance the case of eg = and the symmetric SET, 
Fl = Ftj = F/2. One finds from Eqs. ^ and that 



i{E-2n + ir) 



2E ' A{E -2n + ir/2)^ + - {2SQ)^ 

i{E + 2n + ir) 



A{E + 2Q + iT/2f + F^ - (2 SQ)^ ' 
Upon performing the inverse Laplace transform. 



(25) 



oo+iO 



ailit) 



(26) 



-oo+iO 



and closing the integration contour around the poles of 
the integrand, we obtain for F > 26il and t^l/T 



0"! 



iit)-il/2)oce-^^-^^^^^^'/^sm{2nt). (27) 



Comparing this result with Eq. (fT2)) we find that the 
decoherence rate is 

Fd = 2 (f - v/r2 - 4(5172^ ^-H" (2(5f7)VF . (28) 

For eo 7^ and eo,r <C O the decoherece rate F^ is 
multiplied by an additional factor [1 — (eo/2il)2]. 
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In a general case, ^r, we obtain in the same 

limit {Tf^ fi ^ Sn,) for the decoherence rate: 



Li R 



i + fi^Vi^L+TRr 



(29) 



It is interesting to compare this result with the fluctu- 
ation spectrum of the charge inside the SET, Eq. (jB8|) . 
Appendix B. Wc find 



Td^2{SLORfSQiO), 



(30) 



where ujr = ^ ASt^ + ep is the Rabi frequency. The latter 
represents the energy splitting in the diagonalized qubit 
Hamiltonian. Thus biOR corresponds to the amplitude of 
energy level fluctuations in a single dot. 

Although Eq. pO|) has been obtained for small fluc- 
tuations <5wfl, it might be approximately correct even 
if biOR is of the order of F. It is demonstrated in 
Fig. 4, where we compare a\\(t) and ai^it), obtained 
from Eqs. pTjl and (p2| (solid line) with those from 
Eqs. (|3ap and (H]) (dashed line) for the decoherence rate 
Yd given by Eq. (|30p . The initial conditions correspond to 
CTii(O) = 1 and CTi2(0) = (respectively, Uaa{^) = Yr/T 
and abb{0) =Tr/T). 

In the case of aligned qubit, however. Re cri2{t) — 
Re (712(0), as was explained above. On the other hand, 
one always obtains from ([5a|) and ^ that Re [ai2 {t — > 
00)] = 0. Therefore the phenomenological Bloch equa- 
tions are not applicable for evaluation of Ke[ai2{t)], 
even in the weak coupling limit (besides the case of 
Re [(712(^-0)] -0). 

In the large coupling regime {Sfl 3> F) the phenomeno- 
logical Bloch equations, Eqs. ([5a|) and cannot be 
used, as well. Consider for simplicity the case of e = 
and Tl,r = T/2. Then one finds from Eq. (27) that 
the damping oscillations between the two dots take place 
at two different frequencies, 2n ± yj ((50)^ — (r/2)^, in- 
stead of the one frequency, cjr = 2Q, given the Bloch 
equations. Moreover, Eq. (|5D|) does not reproduce the 
decoherence (damping) rate in this limit. Indeed, one 
obtains from Eq. (28) that the decoherence rate F^ — 2F 
for 6il > F/2, so F^ does not depend on the couphng 
{Sn) at all. 



B. Fluctuation of the energy level 

Consider the SET placed near one of the qubit dots, 
as shown in Fig. 2b. In this case the qubit-SET inter- 
action term is given by Eq. p^ . As a result the energy 
level El will fluctuate under the influence of the fluctua- 
tions of the electron charge inside the SET. The available 
discrete states of the entire system are shown in Fig. 5. 
Using Eqs. (I20p we can write the rate equations, similar 



0.4 
0.2 


-0.2 
-0.4 



Re Oj^it) 



10 20 30 40 50 



(a) 




(b) 



FIG. 4: The occupation probability of the first dot of the qubit 
for e = 2n, Fl = n, Ffl = 2n and Sn = 0.50. The solid line 
is the exact result, whereas the dashed line is obtained from 
the Bloch-type rate equations with the decoherence rate given 
by Eq. dSO}- 





E, + U 





FIG. 5: The available discrete states of the entire system for 
the configuration shown in Fig. 2b. Here U is the repulsion 
energy between the electrons. 



to Eqs. (Ell), 



0"aa = 








— fca), 


(31a) 


CTfcb = 


— r'^CTfeb + 


Y'h'^aa — 


ifloicTbd 


— Cdb), 


(31b) 


O'cc = 


-^Lf^cc + 


^R'^dd — 




- C^ac), 


(31c) 


O'rfd = 


-^R<^dd -1 


- ^Lf^cc - 


'i^oio'db 


- CTbd), 


(31d) 


O'ac = 


— ieofac - 




— CTcc) - 


Tl + F' 
2 





rflF'^cTbd, (31e) 
Tr 



F'„ 



-o-bd 



TLT'i^aac , (31f) 



where F^ are the tunneling rate at the energy Eq + U—. 

Let us assume that F'^^ ^ = F^ j^. Then it follows from 
Eqs. PT|) that the behavior of the charge inside the SET 
is not affected by the qubit, the same as in the previous 
case of the Rabi frequency fluctuations. Also the qubit 
density matrix becomes the mixture ([5]) in the stationary 
state for any values of the qubit and the SET parameters. 
Hence, there is no qubit relaxation in this case either 
(except for the static qubit, fio = 0, and (711(0) 7^ (722(0), 
Eq. dni)). 

Since according to Eq. ([M)) , the probability of finding 
an electron inside the SET in the stationary state is Pi — 
Fi/F, the energy level Ei of the qubit is shifted by PiU . 
Therefore it is useful to define the "renormalized" level 
displacement, e = -\- PiU . 
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As in the previous case we use the Laplace transform, 
a{t) — *■ cr(-E), in order to determine the decoherence rate 
analytically. In the case of = Tji — T/2 and e = we 
obtain from Eqs. ([51]) 



32iE + iT)nl 



(32) 



C/2 - AE(E + iV) 



The position of the pole in the second term of this expres- 
sion determines the decoherence rate. In contrast with 
Eq. ([25|l . however, the exact analytical expression for the 
decoherence rate (Td) is complicated, since it is given by 
a cubic equation. We therefore evaluate in a different 
way, by substituting E — ±2Qq — ij in the second term 
of Eq. ([5^ and then expanding the latter in powers of 
7 by keeping only the first two terms of this expansion. 
The decoherence rate Yd is related to 7 by = 47, as 
follows from Eq. Then we obtain: 



2(r^ + 4f7^) 



for [/< (f72+ri7o)i/2 



In general, if 7^ F^j, one finds from Eqs. ([3T|l that 
Fd = 2U^TlTr/[T{T^ + for U « {Ql + Tno^/^ 

The same as in the previous case, Eq. the deco- 

herence rate in a weak coupling limit is related to the 
fluctuation spectrum of the SET, Sq{oj), Eq. (|B8|l . but 
now taken at a different frequency, uj = 2Qq. The lat- 
ter corresponds to the level splitting of the diagonalized 
qubit's Hamiltonian, un. Thus, 



(34) 



which can be applied also for e 7^ 0. This is illustrated 
by Fig. 6 which shows aii{t) obtained from Eqs. ((3T|) 
and (221) (sohd line) with Eqs. and ^ (dashed hne) 
for the decoherence rate F^ given by Eq. As in 

the previous case, shown in Fig. 4, the initial conditions 
correspond to fTii(O) = 1 and 0-12(0) — (respectively, 
aaa{0) = Fij/F and cr6f,(0) = F/j/F). One finds from 
Fig. 6 that Eq. can be used for an estimation of F^ 
even for t/ ~ F, Hq. 

In contrast with the tunneling-coupling fluctuations, 
Eq. ([HO)) , where the decoherence rate is given by Sq{0), 
the fluctuations of the qubit's energy level generate the 
decoherence rate, determined by the fluctuation spec- 
trum at Rabi frequency, Sq{ujii), Eq. p4p . A similar 
distinction between the decoherence rates generated by 
different components of the fluctuating field, exists in a 
phenomenological description of magnetic resonance^^. 
One can understand this distinction by diagonalizing the 
qubit's Hamiltonian. In this case the Rabi frequency, 
ujR, becomes the level splitting of the qubit's states 
|±) = (|1) ± |2))/a/2 (for e = 0). So in this basis, the 
tunneling-coupling fluctuations correspond to simultane- 
ous fluctuations of the energy levels in the both dots. 



ima.Jt) 



Re O^lt) 




tQo 



(a) 



(b) 



FIG. 6: The probability of finding the electron in the first 
dot of the qubit for e = 2flo, Tl = f^o, Fh = 2^0 and U = 
0.5Slo. The solid line is the exact result, whereas the dashed 
line is obtained from the Bloch-type rate equations with the 
decoherence rate given by Eq. p4)) . 



Since these fluctuations are "in phase" , we could expect 
that the corresponding dephasing rate is determined by 
spectral density at zero frequency. In fact, it looks like 
as fluctuations of a single dot state, considered by Levin- 
son in a weak coupling limit^. On the other hand by 
fluctuating the energy level in one of the dots only, one 
can anticipate that the corresponding dephasing rate is 
determined by the fluctuation spectrum at the Rabi fre- 
quency, ujR, Eq. ([M]) . which is a frequency of the inter-dot 
transitions. 

Since ujji can be controlled by the qubit's levels dis- 
placement, e, the relation (|34p can be implied by using 
qubit for a measurement of the shot-noise spectrum of the 
environmen t^^'^^i"^° . For instance, it can be done by at- 
taching a qubit to reservoirs at different chemical poten- 
tials. The corresponding resonant current which would 
flow through the qubit in this case, can be evaluated via 
a simple analytical expression^ that includes explicitly 
the decoherence rate, Eq. Thus by measuring this 

current for different level displacement of the qubit (eo), 
one can extract the spectral density of the fluctuating 
environment acting on the qubit^. 

Although Eq. for the decoherence rate has been 
obtained by using a particular mechanism for fluctuations 
of the qubit's energy levels, we suggest that this mecha- 
nism is quite general. Indeed, the rate equations pip can 
describe any fluctuating media near a qubit, driven by the 
Boltzmann type of equations. Therefore it is rather nat- 
ural to assume that Eq. p4p would be valid for any type 
of such (classical) environment in weak coupling limit. 
This implies that the decoherence rate is always deter- 
mined via the spectral density of a fluctuating qubit's 
level, whereas the nature of a particular medium inducing 
these fluctuations would be irrelevant. In order to sub- 
stantiate this point it is important to compare Eq. (j34[) 
with the corresponding decoherence rate induced by the 
thermal environment in the framework of the spin-boson 
model. In a weak damping limit this model predictsii^ 
Tj~^ — T^^ = ((7g/2)S'(a;fl) , where go is a couphng of 
the medium with the qubit levels (go corresponds to U in 
our case) and S{ijj) is a spectral density. Using Eq. JIT 
one flnds that this result coincides with Eq. 
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FIG. 7: The probability of finding the electron in the first 
dot of the qubit for e = 0, Fl = Fa = and U, as given 
by Eqs. (|3H) (solid line) and from the Bloch-type equations 
(dashed line) with the decoherence rate given by Eq. (I33|l . 



C. Strong-coupling limit and localization 



Let us consider the limit of fj > {Q^ + Filo)^/^. Our 
rate equation (|3ip are perfectly valid in this region, pro- 
viding only that Eq + U is deeply inside of the potential 
bias, Eq. (dH). We find from Eq. Q that the deco- 
herence rate is not directly related to the spectrum of 
fluctuations in strong coupling limit. In addition, the ef- 
fective frequency of the qubit's Rabi oscillations (w|/^) 
decreases in this limit. Indeed, by using Eqs. ([32]) . ([26|) . 
one finds that the main contribution to f7ii(i), is coming 
from a pole of aii{E), which lies on the imaginary axis. 
This implies that the effective frequency of Rabi oscilla- 
tions strongly decreases when U ^ {Qq + FJIq)^/^. In 
addition, the decoherence rate F^ ^ in the same limit, 
Eq. ([55)1 . As a result, the electron would localize in the 
initial qubit state, Fig. 7. 



The results displayed in this figure show that the solu- 
tion of the Bloch-type rate equations, with the decoher- 
ence rate given by Eq. p3p . represents damped oscilla- 
tions (dashed line). It is very far from the exact result 
(solid line), obtained from Eqs. ([5T|) and corresponding 
to the electron localization in the first dot. The latter is 
a result of an effective decrease of the Rabi frequency for 
large U that slows down electron transitions between the 
dots. Thus such an environment-induced localization is 
different from the Zeno-type effect (unlike an assumption 
of Ref.^^). Indeed, the Zeno effect takes place whenever 
the decoherence rate is much larger then the coupling be- 
tween the qubit's states^i^^. However, the decoherence 
rate in the strong coupling limit is much smaller then the 
coupling r^o . In fact, the localization shown in Fig. 7 
is rather similar to that in the spin-boson modeJii^. It 
shows that in spite of their defferences, both models trace 
the same physics of the back-action of the environment 
(SET) on the qubit. 



BACK-ACTION OF THE QUBIT ON THE 
ENVIRONMENT 



A. Weak back-action effect 

Now we investigate a weak dependence of the width's 
^L,R on the energy U, Fig. 5. We keep only the linear 
term, ^ = TL^ii + aL,RU, by assuming that U is small. 
(A similar model has been considered in^^i^). In con- 
trast with the previous examples, where the widths have 
not been dependent on the energy, the qubit's oscillation 
would affect the SET current and its charge correlator. 
A more interesting case corresponds to 7^ an. Let us 
take for simplicity q;l = and an = a 0. 

Similarly to the previous case we introduce the "renor- 
malized" level displacement, e ^ eq — {Tl/T)U, where 
e = corresponds to the aligned qubit. Solving Eqs. ([5T|) 
in the steady-state limit, a = a{t oo), and keeping 
only the first term in expansion in powers of J7, we find 
for the reduced density matrix of the qubit, Eqs. 



ca 



2 4Ffl 

ano(l + caU) 
2f^ 



2Tr 
1 ae 

- A 

2 



4F 



) 



(35) 



where c = (ae - 2F)/(4FflF). It follows from Eqs. JSH) 
that the qubit's density matrix in the steady-state is no 
longer a mixture, Eq. ([5]) . Indeed, the probability to 
occupy the lowest level is always larger than 1/2 and 
CT12 7^ 0. This implies that relaxation takes place to- 
gether with decoherence. The ratio of the relaxation and 
decoherence rates is given by the off-diagonal terms of 
the reduced density matrix of the qubit. For e = one 
finds from Eq. ^ that Fd/F^ = a^^ " 2. 

In order to find a relation between the decoherence 
and relaxation rates, F^.r, and the fluctuation spectrum 
of the qubit energy level, Sq{lo), we first evaluate the 
total damping rate of the qubit's oscillations (7). Using 
Eq. (fT^ we find that this quantity is related to the deco- 
herence and relaxation rates by 7 = (F^ -I- 2Fr)/4. The 
same as in the previous case the rate 7 is determined by 
poles of Laplace transformed density matrix a{t) — > (j{E) 
in the complex E'-plane. Consider for simplicity the case 
of e = and F^ = Vr = F/2. Performing the Laplace 
transform of Eqs. pT|) we look for the poles of aii{E) at 
E — ±2ilo — *7 by assuming that 7 is small. We obtain 



2(F2 + Afll) 



4nl 



2 (F2 + 4112) 



(36) 



for U < V,o. 

Now we evaluate the correlator of the charge inside the 
SET, Sq{lu) which induces the energy-level fluctuations 
of the qubit. Using Eqs. ((3T|) and (|B6|) we find. 



F 



2(F2 +a;2) 



aU- 



F2 



4(F2+w2)^ 



(37) 
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for all < r. Therefore in the hmit of C7 < fig and a ?7 < 
r the total damping rate of the qubit's oscillations is 
directly related to the spectral density of the fluctuations 
spectrum taken at the Rabi frequency, 



1^. 



(38) 



This represents a generalization of Eq. (|34p for the case 
of a weak back-action of qubit oscillations on the spectral 
density of the environment. As a result, the qubit dis- 
plays relaxation together with decoherence. It is remark- 
able that the total qubit's damping rate is still given by 
the fluctuation spectrum of the SET (environment) mod- 
ulated by the qubit. Note that Eq. ([55)) can be apphed 
only if the modulation of the tunneling rate through the 
SET (tunneling current) is small a ?7 <C F, in addition to 
a weak distortion of the qubit {U ^ fto)- 

In the case of strong back-action of the qubit on the 
environment the decorerence and relaxation rates of the 
qubit are not directly related to the fluctuation spectrum 
of the environment, even if the distortion of the qubit is 
small. This point is illustrated by the following example. 



B. Strong back-action 

Until now we considered the case where Eq + U <ti 
jjLL , so that the interacting electron of the SET remains 
deeply inside the voltage bias. If however, the interaction 
U between the qubit and the SET is such that Eq + 
U ^ fiL, the qubit's oscillation would strongly affect the 
fluctuation of charge inside the SET. Indeed, the current 
through the SET is blocked whenever the level Ei of the 
qubit is occupied. Fig. 8. In fact, this case can be treated 
with small modiflcation of the rate equations (jSip . if only 
//i — i?o > F and Eq + U — fj-L T, where Eq is a level 
of the SET carrying the current. 

The corresponding quantum rate equations describing 
the system are obtained directly from Eqs. ([20]) . As- 
suming that the widths r^,/;, are energy independent we 
findi^ 

CTaa = (Fl + TR)abb - ifto{(Tac - (Tea), (39a) 

o-bb = -(Fij + TL)(Jbb - i^oio-bd - (Jdb), (39b) 

O-cc = -FlCTcc + Tfiadd ~ i^oio-ca ~ Oac), (39c) 
Odd = —^R(Tdd + FlCTcc — iVLQ{(Jdb — CTbd), (39d) 

Fl 



CTa 



ifio(o'aa 



-On 



+ ^RPbd, (39e) 
Obd = -i(eo + U)obd - i^oiobb ~ Odd) 

Fl 



rR + ^]obd. (39f) 



Solving Eqs. (|39|) in the stationary limit, a = a{t —i- 
oo) and introducing the "renormalized" level displace- 
ment, e = Co — ?7Fl/(2F), we obtain for the qubit's den- 



□ 






FIG. 8: The available discrete states of the entire system 
when the electron-electron repulsive interaction U breaks off 
the current through the SET. 



sity matrix, Eqs. in the steady state: 



Oil 



Ol2 



1 8eU 

i2uno 

16e2 -I- 8C/e -I- ASQ^ + 9{U^ + T^) ' 



,(40a) 
(40b) 



where for simplicity we considered the symmetric case, 
Fl = Fi^ = F/2. It follows from Eqs. (gOl) that similarly 
to the previous example, the qubit's density matrix is no 
longer a mixture ^ . The relaxation takes place together 
with decoherence in this case too. 

Let us consider weak distortion of the qubit by the 
SET, f/ < r^Q. Although the values of U are restricted 
from below {U 3> T + ^l — Eq), this limit can be achieved 
if the level Eq is close to the Fermi energy, providing only 
that /iL — -Eo 3> F, and F <C C/. Now we evaluate oii{t) 
with the rate equations (j39p and then compare it with 
the same quantity obtained from the Bloch equations, 
Eq. in]), where F^^^ are given by Eqs. jSUand The 
corresponding charge-correlator, Sq{u!ii), is evaluated by 
Eqs. (jB6[) and ([55)) . As an example, we take symmetric 
qubit with aligned levels, e = 0, Fl = Fj^ = O.OSfio and 
U = O.Sr^o- The decoherence and relaxation rates, corre- 
sponding to these parameters are respectively: Td/^la — 
0.0038 and Fjl^o 0.00059. 

The results are presented in Fig. 9a. The solid line 
shows aii{t), obtained from the rate equations (|39p . 
where the dashed line is the same quantity obtained from 
Eq. ini). We find that Eq. ^ (or ([38])) underestimates 
the actual damping rate of oii{t) by an order of mag- 
nitude). This lies in a sharp contrast with the previous 
case, where the energy level of the SET is not distorted 
by the qubit, F'^ ^ = FL,i?, Fig. 5. Indeed, in this case 
crii(i) obtained Eq. ([12]) with F^ given by Eq. ([34]) and 
Tr = 0, agrees very well with that obtained from the rate 
equations ([5T]) . as shown in Fig. 9b. 

Such an example clearly illustrates that the decoher- 
ence is not related to the fluctuation spectrum of the 
environment, whenever the environment is strongly af- 
fected by the qubit, even if the coupling with a qubit 
is small. This is a typical case of measurement, corre- 
sponding to a noticeable response of the environment to 
the qubit's state (a "signal"). 
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(b) 

FIG. 9: (a) The probability of finding the electron in the first 
dot of the qubit for e = 0, Fl ^Tr = O.OSrio and U = 0.5flo. 
The solid line is obtained from Eqs. (|39|) . whereas the dashed 
line corresponds to the Eq. (|12|l with Td given by Eq. (|34}; 
(b) the same for the case, shown in Fig. 5, where the solid 
line corresponds to Eqs. pT|l. 



VI. SUMMARY 

In this paper we propose a simple model describing 
a qubit interacting with fluctuating environment. The 
latter is represented by a single electron transistor (SET) 
in close proximity of the qubit. Then the fluctuations 
of the charge inside the SET generate fluctuating field 
acting on the qubit. In the limit of large bias voltage, the 
Schrodinger equation for the entire system is reduced to 
the Bloch-type rate equations. The resulting equations 
are very simple, so that one can easily analyze the limits 
of weak and strong coupling of the qubit with the SET. 

We considered separately two different cases: (a) there 
is no back-action of the qubit on the SET behavior, so 
that the latter represents a "pure environment" ; and (b) 
the SET behavior depends on the qubit's state. In the 
latter case the SET can "measure" the qubit. The setup 
corresponding to the "pure environment" is realized when 
the energy level of the SET carrying the current lies 
deeply inside the potential bias. The second (measure- 
ment) regime of the SET is realized when the tunnel- 
ing widths of the SET are energy dependent, or when 
the energy level of the SET carrying the current is close 
enough to the Fermi level of the corresponding reservoir. 
Then the electron-electron interaction between the qubit 
and the SET modulates the electron current through the 
SET. 

In the case of the "pure environment" ("no- 
measurement" regime) we investigate separately two dif- 
ferent configurations of the qubit with respect to the 
SET. In the first one the SET produces fluctuations of 
the off-diagonal coupling (Rabi frequency) between two 
qubit's states. In the second configuration the SET pro- 
duces fiuctuations of the qubit's energy levels. In the 
both cases we find no relaxation of the qubit, despite 
the energy transfer between the qubit and the SET can 
take place. As a result the qubit always turns asymptot- 
ically to the statistical mixture. We also found that in 
both cases the decoherence rate of the qubit in the weak 
coupling limit is given by the spectral density of the cor- 



responding fluctuating parameter. The difference is that 
in the case of the off-diagonal coupling fluctuations the 
spectral density is taken at zero frequency, whereas in 
the case of the energy level fluctuations it is taken at the 
Rabi- frequency. 

In the case of the strong coupling limit, however, the 
decoherence rate is not related to the fluctuation spec- 
trum. Moreover we found that the electron in the qubit 
is localized in this limit due to an effective decrease of 
the off-diagonal coupling. This phenomenon may resem- 
ble the localization in the spin-boson model in the strong 
coupling limit. 

If the charge correlator and the total SET current are 
affected by the qubit (back-action effect), we found that 
the off-diagonal density-matrix elements of the qubit sur- 
vive in the steady-state limit and therefore the relax- 
ation rate is not zero. We concentrated on the case of 
weak coupling, when the Coulomb repulsion between the 
qubit and the SET is smaller then the Rabi frequency. 
The back-action of the qubit on the SET, however, can 
be weak or strong. In the first case we found that the 
total damping rate of the qubit due to decoherence and 
relaxation is again given by the spectral density of the 
SET charge fluctuations, modulated by the qubit. This 
relation, however, is not working if the back-action is 
strong. Indeed, we found that the damping rate of the 
qubit in this case is larger by an order of magnitude than 
that given by the spectral density of the corresponding 
fluctuating parameter. 

This looks like that in the strong back-action of the 
qubit on the SET the major component of decoherence is 
not coming from the fluctuation spectrum of the qubit's 
parameters only, but also from the measurement "sig- 
nal" of the SET. On the flrst sight it could agree with 
an analysis of Ref.— , suggesting that the decoherence 
rate contains two components, generated by a measure- 
ment and by a "pure environment" (environmental fluc- 
tuations). The latter therefore represents an unavoid- 
able decoherence, generated by any environment. Yet, 
in a weak coupling regime such a separation seems not 
working. In this case the damping (decoherence) rate is 
totally determined by the environment fluctuations, even 
so modulated by the qubit. 

Although our model deals with a particular setup, it 
bears the main physics of a fluctuating environment, act- 
ing on a qubit. Indeed, the Bloch-type rate equations, 
which we used in our analysis have a pronounced phys- 
ical meaning: they relate the variation of qubit param- 
eters with a nearby fluctuating field described by rate 
equations. A particular mechanism, generated this field 
should not be relevant for an evaluations of the deco- 
herence and relaxation rates, but only its fluctuation 
spectrum. Indeed, in the weak coupling limit our re- 
sult for the decorence rate coincides with that obtained 
in a framework of the spin-boson model. Thus our model 
can be considered as a generic one. Its main advantage is 
that it can be easily extended to multiple coupled qubits. 
Such an analysis would allow to determine how decoher- 
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ence scales with number of qubits^, which is extremely 
important for a realization of quantum computations. 

In addition, our model can be extended to a more 
complicated fluctuating environments, such as containing 
characteristic frequencies in its spectrum. It would for- 
mally correspond to a replacement of the SET in Fig. 2 
by a double-dot (DD) coupled to the reservoirs^. All 
these situations, however, must be a subject of a sepa- 
rate investigation. 
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FIG. 10: Resonant tunneling through a single dot. ^l,r are 
the Fermi energies in the collector and emitter, respectively. 



amplitudes h{E): 

Eh{E) -Y,^ME) ^ I (A2a) 

I 

(E + El ^ Eo)bMiE) ~ riME) 

-J2^rbiriE)^0 (A2b) 

r 

{E + El - Er)blr{E) ^ nrbol{E) 

-Y,^i'iowriE) = (A2c) 
{E + Ei+ El, -Eo- Er)boii'r{E) - ni,bi,{E) 

+ flihriE) - J2 ^r~blVrr'{E) = (A2d) 



APPENDIX A: QUANTUM-MECHANICAL 
DERIVATION OF RATE EQUATIONS FOR 
QUANTUM TRANSPORT 



Consider the resonant tunneling through the SET, 
shown schematically in Fig. 10. The entire system is 
described by the Hamiltonian i?sET, given by Eq. (fH)) . 
The wave function can be written in the same way as 
Eq. (fT7)) . where the variables related to the qubit are 
omitted, 



*(t)) = \b{t) + ^ 6o/(t)4c; + ^ri{t)cl 



Cl 



^ borU'{t)cl 



l,r 



t t 
Cic/Q' 



l<l',r 



|0). (Al) 



Substituting \^{t)) into the time-dependent Schrodinger 
equation, idt\^{t)) — iJgExl and performing the 
Laplace transform, b{E) = exp{iEt) b{t)dt, we obtain 
the following infinite set of algebraic equations for the 



(The r.h.s of Eq. (jA2aP reflects the initial condition.) 

Let us replace the amplitude b in the term ^ fib of each 
of the equations (IA2|) by its expression obtained from the 
subsequent equation. For example, substituting boi{E) 
from Eq. (jA2bjl into Eq. (|A2aP we obtain 



Ei-Eo 



E 

l,r 



b{E) 



E + Ei-Eo 



bir{E) 



(A3) 



Since the states in the reservoirs are very dense (contin- 
uum), one can replace the sums over I and r by integrals, 
for instance J pL{Ei)dEi , where pl{Ei) is the 

density of states in the emitter, and £7;^^ — *■ ^L.R{Ei^r)- 
Consider the first term 



^1 = 



^l{Ei) 

E + Ei-Eo 



PL{Ei)dEi 



(A4) 



where A is the cut-off parameter. Assuming weak en- 
ergy dependence of the couplings 0.lm and the density of 
states pl,r, we find in the limit of high bias, /x^ = A — > c» 

Si^~ninl{Eo~E)pL{EQ-E) = -i^. (A5) 
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Consider now the second sum in Eq. (|A3|) 



S2 



PB.{Er)dEr 



-A 



X / pL[Ei)dEi, (A6j 



A 



E + Ei- Ei 



^0 



where we replaced bir{E) by b{E, Ei, E^) and took /i^ = 
A, pn = —A. In contrast with the first term of Eq. (IA3p . 
the ampUtude 6 is not factorized out the integral (|A6ll . 
We refer to this type of terms as "cross-terms" . Fortu- 
nately, all "cross-terms" vanish in the limit of large bias, 
A — !■ CX3. This greatly simplifies the problem and is very 
crucial for a transformation of the Schrodinger to the rate 



J 



equations. The reason is that the poles of the integrand 
in the £^;(i?j.)-variable in the "cross-terms" are on the 
same side of the integration contour. One can find it by 
using a perturbation series the amplitudes b in powers of 
fi. For instance, from iterations of Eqs. (jA2[) one finds 



biE,Ei,Er 



EiE + El - Er)iE + El - Eo) 



(A7) 



The higher order powers of have the same structure. 
Since E ^ E + ie in the Laplace transform, all poles of 
the amplitude b{E, Ei, Er) in the i?;-variable are below 
the real axis. In this case, substituting Eq. (jA7p into 
Eq. ((X6|) we find 



lim 

A^oo 



-A 



{E + ie){E + E0-E1+ ieY{E + Eq - Er + ie) 

I 



dEi^Q, 



(A8) 



Thus, 5*2 in the limit of /i^ 00, /i/j — s- —00. 

Applying analogous considerations to the other equa- 
tions of the system (|A2[) , we finally arrive at the following 
set of equations: 

{E + iTL/2)b{E) =i (A9a) 
{E + Ei-Ea+iTR/2)bm{E) 

- VLib{E) ^ (A9b) 

{E + Ei^Er+lTLl2)bir{E) 

-VLrbQi{E)^Q (A9c) 
{E + Ei+ Ev -Eo- Er + iTR/2)bm'r{E) 

-ni.bir{E)+nibvr{E)={) (A9d) 



l...,r... •' 



l...,r.. 



Eqs. (jA9[) can be transformed directly to the reduced 

density matrix crj"/'" \t), where j = 0, 1 denote the state 
of the SET with an unoccupied or occupied dot and n de- 
notes the number of electrons which have arrived at the 
collector by time t. In fact, as follows from our derivation, 
the diagonal density-matrix elements, j — j'and n = n' , 
form a closed system in the case of resonant tunneling 
through one level. Fig. 10. The off-diagonal elements, 
j 7^ j' , appear in the equation of motion whenever more 
than one discrete level of the system carry the transport 
(see Eq. (|20l) . Therefore we concentrate below on the di- 



agonal density-matrix elements only, aQ^\t) = cT^'"''{t) 
and a[^\t) = crj"'"''(t). Applying the inverse Laplace 



(n,n) / 



transform on finds 



dEdE' ~ ~ fp'\J{E'-E)t 
—bi,,,^,,,{E)bi_^^_iE)e'- 




kl---r--XE)b;i,,,^,,,iE')e^(^'-^^' 



(AlOa) 

(AlOb) 



Consider, for instance, the term <J^ii (t) — J2i l&o/(*)P- E ^ E' we obtain 
Multiplying Eq. (jA9bp by b^i{E') and then subtracting 
the complex conjugated equation with the interchange 



dEdE' 
47r2 

-21mY,^iboiiE)b*{E' 



J2[iE' - E - iTRyboiiE)b*i{E') 
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Using Eq. (jAlObp one easily finds that the first inte- Here SQ{t) = CQ{t)co{t) — q and q = Pi = Pi(t ^ oo) is 



gral in Eq. (jAlip equals to 
substituting 



a['l>{t) + rna\\>{t)]. Next, 



(0), 



E + E1-E0 + irR/2 



(A12) 



from Eq. (jA9bp into the second term of Eq. (|A11[) . and 
replacing a sum by an integral, one can perform the Ei- 
integration in the large bias limit, /Lt^ ^ 00, /Lt/j — > — 00. 
Then using again Eq. (jAlObp one reduces the second 

term of Eq. (Ku\ to irLcr^o^(t). Finally, Eq. (|XTT|) reads 



The same algebra can be applied for all other am- 
plitudes ba{t). For instance, by using Eq. (jAlOap one 
easily finds that Eq. (|A9c[) is converted to the following 
rate equation a^ol\t) = -ri,CT^o^(t) + Tna^^^'it). With 
respect to the states involving more than one electron 
(hole) in the reservoirs (the amplitudes like hm'r{E) and 
so on) , the corresponding equations contain the Pauli ex- 
change terms. By converting these equations into those 
for the density matrix using our procedure, one finds the 
"cross terms", like J2^i^i'r{E)fli'b'^^{E'), generated by 
Eq. (jA9d|l . Yet, these terms vanish after an integration 
over Ei(r ) in the large bias limit, as the second term in 
Eq. (|A3p . The rest of the algebra remains the same, 
as described above. Finally we arrive at the following 
infinite system of the chain equations for the diagonal 



elements, ctqq'' and a\'i , of the density matrix, 



An) 



(0--ri4°)(t), 



^(0) 



(0), 



(Al3a) 
(A13b) 
(Al3c) 
(A13d) 



Summing over n in Eqs. (jA13|) we find for the reduced 



density matrix of the SET, a{t) 
lowing "classical" rate equations, 



E„CT(")(t), the fol- 



^Too(i) = ~TLaooit)+TRaii{t) (A14a) 
CTii(i) = TLaoo{t)~rRaii{t) (A14b) 

These equations represent a particular case of our general 
quantum rate equations (|20[) , which are derived using the 
above described technique^i^^. 



APPENDIX B: CORRELATOR OF ELECTRIC 
CHARGE INSIDE THE SET. 

The charge correlator inside the SET is given by 
Sq{u)) = Sq{uj) + Sq{-uj), where 



'5'q(w) 



{5Q{0)6Qit))e"^'dt. 



(Bl) 



the average charge inside the dot. Since the initial state, 
t = in Eq. (jBip corresponds to the steady state, one 
can represent the time-correlator as 

{6Q{0)SQ{t)) = P,{0){q - q){{Q,{t)) - q) , (B2) 

9=0,1 

where Pq{0) is the probability of finding the charge q = 
0, 1 inside the quantum dot in the steady state, such that 
Pi(0) = q and Po(0) = 1-9, and {Qg{t)) = P^'^t) is 
the average charge in the dot at time t, starting with the 



initial condition Pj''' (0) = q. 
Eq. (jBip we finally obtain 



Substituting Eq. (|B2]) into 



SQ{co)=q{l-q)[P^'\uj) 



(B3) 



where P^'^\uj) is a Laplace transform of P['^\t). These 
quantities are obtained directly from the rate equations, 
such that q — Ubb + &dd and p['^\ijj) — all\uj) + 
a'fJ^Lo), where a ~ a{t 00) and (T^^^(a;) is the 
Laplace transform cr^'?-' (t) with the initial conditions 
corresponding to the occupied (q = 1) or unoccupied 
(q — 0) SET. In order to find these quantities it is use- 
ful to rewrite the rate equations in the matrix form, 
&{t) — M(j{t), representing (j{t) as the eight- vector, 

= {(^aai(^bb,(^cc,<yddT'^ac,<yca,(^bd,(^db} and M aS the 

corresponding 8 x 8-matrix. Applying the Laplace trans- 
form we find the following matrix equation. 



(B4) 



where / is the unit matrix and ct'-'^-' (0) is the initial con- 
dition for the density-matrix obtained by projecting the 
total wave function (|17p on occupied (g = 1) and unoc- 
cupied (q = 0) states of the SET in the limit oi t 00, 

ct(i)(0) = AAi{0,CT6b,0,o-dd,0,0,CTM,o-d6}, (B5a) 

aW(0) = AAo{(7aa,0,acc,0,aac,^Tca,0,0}, (B5b) 

and A/i — 1/q and A/q — 1/(1 — g) are the corresponding 
normalization factors. Finally one obtains: 

Sq{lu) = 2q{l - q)Re [aH^co) + a<;lJ (co) 

-^il\^)~^'dJi^)]- (B6) 

In the case shown in Fig. 2 one finds from Eqs. (|2ip or 
Eqs. (EID for T'^^R = Tl,r that aac = (Tm = 0, q = Fl/F 

and a^^^(aj) -I- ^'f]{uj) — Pi''\uj). The latter equation is 
given by 



{iu;^T)P^'\lo) = -q+ — 



Substituting Eq. (|B7)l into Eq. (|B3|l one obtains: 

2TlTr 



F(cj2-KF2) 



(B7) 



(B8) 
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Obviously, for a more general case when ^ ^l,r, 
or when the electron-electron interaction excites the elec- 
tron inside the SET above the Fermi level, Fig. 8, the ex- 



pressions for Sq{uj), obtained from Eq. (|B6|) have a more 
complicated than Eq. (|B8[) . 
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